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Abstract 

This lecture is a sketch of the physics of the cosmic microwave background. 
The observed anisotropy can be divided into four main contributions: varia- 
tions in the temperature and gravitational potential of the primordial plasma, 
Doppler effect from its motion, and a net red/blueshift the photons accumulate 
from traveling through evolving gravitational potentials on their way from the 
primordial plasma to here. These variations are due to primordial perturba- 
tions, probably caused by quantum fluctuations in the very early universe. The 
ongoing Planck satellite mission to observe the cosmic microwave background 
is also described. 

1 Introduction 

The cosmic microwave background (CMB) is radiation that comes from the early universe. In the early 
universe, ordinary matter was in the form of hot hydrogen and helium plasma which was almost homoge- 
neously distributed in space. Almost all electrons were free. Because of scattering from these electrons, 
the mean free path of photons was short compared to cosmological distance scales: the universe was 
opaque. As the universe expanded, this plasma cooled, and first the helium ions, then also the hydrogen 
ions captured the free electrons: the plasma was converted into gas and the universe became transparent. 
After that the photons of the thermal radiation of this primordial plasma have travelled through the uni- 
verse and we observe them today as the CMB. The CMB is close to isotropic, i.e., the microwave sky 
appears almost equally bright in every direction. With sensitive instruments we can, however, see small 
variations, the CMB anisotropy. 




Fig. 1: The sky at optical and microwave wavelengths: A sky map of the first two week of observations by the 
Planck satellite at the 70 GHz frequency, superimposed on an optical image of the sky. From Ref. Credit: 
ESA, LFI & HFI Consortia (Planck), Optical image: Axel Mellinger 



There is also "foreground" microwave radiation that comes from astiophysical sources, our own 
galaxy and other galaxies. In Fig. [T] we see the radiation from the Milky Way as a horizontal red band in 
the microwave image, whereas further away from the galactic plane we see variations in the intensity of 
the CMB. The foreground can be separated from the CMB by measuring at several frequencies, since it 
has a different electromagnetic spectrum. 

The formation of helium and hydrogen atoms is called recombination, although in this context it 
is a misnomer, since this was the first time the ions and electrons formed atoms. The related increase 
of the photon free mean path beyond cosmological distance scales is called photon decoupling. This 
happened when the age of the universe was about 380 000 years old. At this time there were small 
density variations, about one part in ten thousand in the primordial plasma/gas. After photon decoupling, 
the over-densities began to grow by gravitational attraction and eventually led to the formation of galaxies 
and stars hundreds of millions of years later. 




Fig. 2: Temperature anisotropy of the CMB according to 7 years of measurements by the WMAP satellite. This 
is a false-color image, where yellow and red indicate hotter than average, and blue colder than average. From 
Ref. El. Credit: NASA / WMAP Science Team. 

When looking at the CMB we are thus looking at the 380 000 year old early universe. We see 
those distant parts of the universe from which it has taken the whole remaining part of the history of the 
universe for the light to travel from there to here. The observed small variations in the CMB reflect the 
small density variations at that time. See Fig.[2l 

Because of the finite speed of light, everything we see lies on our past light cone (see Figure O. 
The intersection of our past light cone with the time of photon decoupling forms a sphere, which we 
call the sphere of last scattering. It is this sphere that we observe when we observe the CMB: we 
see each photon coming from the location where it last scattered from an electron. When the photons 
travel from the last scattering sphere to here they are redshifted by the expansion of the universe. The 
universe has expanded by a factor of 1100 since last scattering, and therefore the photon wavelengths 
have been stretched by that factor. Photons decoupled when the temperature of the primordial plasma/gas 
was about 3000 K, and therefore the photons had then a blackbody spectrum with that temperature. 
When all wavelengths of a blackbody spectrum are stretched by the same factor, the spectrum remains 
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blackbody, but its temperature falls with the same factor. The observed mean temperature of the CMB is 
To = 2.725 ±0.001 K today |3|. 

However, because of the inhomogeneity of the universe, photons coming from different directions 
have suffered slightly different redshifts, which is another contribution to the observed CMB anisotropy. 
Thus the variation dT{6, (j)) of the observed temperature T{9, (j)) = 2.72bK + 5T{9, (/)) can be divided 
into two contributions, JTintr that is due to inhomogeneous conditions at the last scattering surface, and 
(5Tjour that arises as the photons travel from the last scattering sphere to here. 




Fig. 3: A spacetime diagram of our past light cone. 

An important thing of the anisotropy ST{9, 4>) is that it is small. The root-mean-square variation 
is about lOO/uK, or 

ST 

— -4x10-^ (1) 
Jo 

While this makes observing this anisotropy very difficult, it simplifies understanding and calculating the 
physics that causes this anisotropy: The primordial density perturbations 6p that are responsible for this 
anisotropy must have also been very small, and we can calculate their evolution using first-order perturba- 
tion theory around a homogeneous and isotropic model of the universe, the so-called background model. 
The deviations from this background model are small, so we can ignore any products of two or more 
such small quantities. This makes the equations linear, so that they can be easily Fourier transformed, 
and lead to a direct relation between initial and final values. 

CMB was discovered by Penzias and Wilson IH in 1964, using a microwave antenna at Bell 
Laboratories in Holmdel, New Jersey. The CMB anisotropy was first measured by the COBE satellite [51 
in 1992, and much more accurate measurements have later been taken by the WMAP satellite and are 
currently being taken by the Planck satellite. 

In this lecture I sketch our present understanding how the CMB anisotropy arises. The relevant 
physics involves quantum field theory in curved spacetime (for the generation of primordial perturbation) 
and general relativistic perturbation theory (for their evolution and effect on the CMB), and it is not 
possible to give a full account in this short lecture. However, many parts of the relevant physics are 
relatively easy to understand, and I try to present those here; for the other parts I just have to give results 
without derivation, in an attempt to present a continuous story. I also give a short overview of the ongoing 
Planck satellite mission to observe the CMB. 
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2 Background Universe 

The background model is the Friedmann-Robertson- Walker (FRW) universe, where the spacetime metric 
is 

ds^ = -dt^ + a{tf [dx^ + dy^ + dz^) (2) 

(we assume here the flat FRW universe, since it agrees well with observations, and simplifies the discus- 
sion). Here a = a{t) is the scale factor that describes the expansion of the universe. The coordinates x, 
y, z are comoving coordinates that are related to actual distances via the scale factor. The expansion rate 
is given by the Hubble parameter 

mt) = i| (3) 

a at 

whose present value Hq = H{to) is called the Hubble constant (Iq denotes the present time). The 
inverse of the Hubble parameter is called the Hubble distance. It gives the distance over which 
causal effects can act in a cosmological time scale; in cosmology it is also called the horizon scale. 
Today it is comparable to the size of the observable universe, but at early times it was much smaller than 
how far we can see today. 

We denote the homogeneous quantities of the background model with an overbar, e.g., p = p{t) 
is the energy density, p = p{t) is the pressure, and so on. The evolution of the background universe is 
governed by the Friedmann equations 



- '-fp (4) 



1 d^a AttG 



ip+m (5) 



a dfi 

where G is the gravitational constant. 

There are at least four different kinds of energy and/or matter, whose energy density makes an 
important contribution to the universe: photons (the CMB), neutrinos, baryonic matter, and cold dark 
matter: 

P = P-i^ Pu + Pb + Pcdm = X] 

P = Pj+Pu+Pb+Pcdm= P'r/^ + Pu/3+Pb = '^Pi (7) 

(Baryonic matter refers to ordinary matter made out of nucleons and electrons.) The early universe was 
radiation-dominated, most of the energy was in the form of photons, neutrinos (and even earlier, also 
other relativistic particles), so that p ~ ^p, leading to an expansion law a oc t^/^. Later the universe 
became matter-dominated, so that^ <C p, leading to an expansion law a oc t^/^. The transition from 
radiation domination to matter domination happened before photon decoupling. 

It appears that a few billion years ago the expansion of the universe began to accelerate, imply- 
ing that a fifth energy component, called dark energy, with negative pressure, p < —^p, had become 
dominant. 



3 The Perturbed Universe 

In the FRW universe there are no stars or galaxies, and no anisotropy of the CMB, since it is completely 
homogeneous and isotropic. The early part of the history of the universe, when deviations from homo- 
geneity were small, and also large scales at later times, can be discussed using perturbation theory. The 
metric of this "perturbed" universe can be written as 

ds'^ = -(1 + 2^)dt^ + 0^(1 - 2*)(dx2 + dy'^ + dz^) (8) 
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where ^{t, x,y, z) may be called the gravitational potential, since in the Newtonian limit of general 
relativity, it indeed becomes the Newtonian gravitational potential due to density perturbations: an over- 
density causes a negative an under-density a positive ^. 

The energy densities and pressures of the different components have perturbations, 

Pi = {1 + 6i)pi and = + Spi (9) 

Since the background model was isotropic, there was no net flow of energy, but in the perturbed universe 
we have also a mean velocity vt of each energy component with respect to the coordinate system. The 
"fluid" perturbation variables 6i, 6pi, and Vi give a sufficient description of the perturbations in each 
component for as long as they can be approximated as perfect fluid, i.e., for as long as the mean free 
paths of their particles are short compared to the distance scales we are interested in. In the perfect fluid 
approximation, the two metric perturbations become equal, <I> = 'f. 

After they have decoupled, the fluid description is not sufficient to describe the evolution of neu- 
trinos and photons. Neutrinos decouple early, during the first second of the history of the universe. After 
that, the neutrino contribution causes a ~ 10% difference between $ and until the universe becomes 
matter-dominated. 



4 Photon Distribution Function 

For the discussion of CMB anisotropy, we need a more detailed description of the photons, given by the 
photon distribution function f(t,x,q), defined so that at a given time t there are 

dN=-^f{t,x,q)dVd^q (10) 

photons with momentum in a range d^q around the value q within a volume dV around location x. (The 
factor 2 comes from the two polarization states of photons.) We divide the photon momentum q into its 
magnitude q (photon energy) and direction q (a unit vector), q = qq. 

In the background model, photons have the blackbody spectrum 

^~^'^^^ = ^Ut^, (11) 

where T{t) is the homogeneous temperature of the distribution. 
In the perturbed universe we write 

f = f + Sf^ i ^ (12) 

defining a momentum-dependent relative temperature perturbation @{t,x,q). Any function f{t,x,q) 
can be written in this form, but the important point is that to 1st order in perturbation theory, Q does not 
develop any dependence on photon energy q. 

Thus we have a radiation temperature perturbation called the brightness function 

e = e{t,x,q) (13) 

which depends just on time, location, and photon direction. Indeed, Fig.|2]is a plot of the measured CMB 
brightness function at our location and time as a function of the observation direction —q. However, to 
predict its properties from theory, we need to follow how @{t,x, q) evolves from early times. 

From Eqs. (fTTT) and (fT2]) the relation between 5f and © is 

Sf = -q%Q (14) 
oq 
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We can integrate the brightness function over the direction angles to get various momenta, or 
multipoles, of the photon distribution. The three lowest ones are called the local monopole, dipole, and 
quadrupole of the photon perturbation: 

Qo{t,x) = ^ I Q{t,x,q)dn 
4vr J 

Qi{t,x) = ^ I qQ{t,x,q)dn 
4vr J 

Q^{t,x) = 1. 1 (^q^q^ -^S,j^eit,x,q)dn (15) 

(we denote the components of x and qby x* = {x^,x'^,x^) and q^ = {q^, q^ , q^)). The monopole and 
dipole of the photon distribution give directly the photon density and velocity perturbation: 

6^ = 490 and = SOi (16) 



5 Boltzmann Equation 

From statistical physics we have the Liouville theorem: If there are now collisions between the particles, 
their distribution function will be constant along any particle trajectory in phase space: 





dt 



dl 
dt 



df dx" 
'dx^ltt 



d£d(f_ 
'dq^~dt 



(17) 



(we sum over repeated indices), where the total derivatives dx'^/dt and dq^ /dt refer to the motion of 
the particle. Collisions between particles modify the equation by adding a collision term C[f], which 
depends on the distribution function: 

| = C[/1. (18) 

Actually, it depends on the distribution functions of all species of particles the photons may collide with, 
in addition to the photon distribution function. 

In curved spacetime, photons travel on lightlike geodesies. The photon coordinate velocity is given 

by 

dx^ _ q" 
dt a 

The evolution of photon momentum and energy is obtained from the geodesic equation of general rela- 
tivity, which can be derived from the metric. For the photon energy it gives 



(19) 



dq 

Tt^'' 

To first order in perturbation theory 



-H 



df dq' 
dq^~dt 



qi 

a dx'- dt 



dldq 
dq dt 



(20) 



(21) 



i.e., we care only about the dependence of the distribution function on, and the change of, photon energy, 
not photon direction. This is because in the background model, the distribution function (fTTl) does not 
depend on the photon direction, and the photons do not change direction (except in collisions). Thus 
both factors on the left-hand side of Eq. (|2T1) have only a small "first-order" direction dependence, so that 
the direction dependence of the product is "second-order small". Thus the photon Boltzmann equation 
becomes 



4f 
dt 



a dx^ dq 



dt 



q' 5$ 
a dx^ 



+ 



d-^ 
'dt 



C[f] 



(22) 
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Of the five terms in Eq. (122] ). the first two aie just kinematics: the distribution function changes 
in time, since photons move in and out of volume elements due to their velocity. The third term gives 
the change in photon energy (redshift) due to the overall expansion of the universe. The two last terms 
give the effect of spacetime perturbations: the gravitational redshift due to a gradient in gravitational 
potential, and the effect of local variations in the expansion rate. 

We can separate Eq. (l22l) into a background equation 







(23) 



dt dt dq 

(the effect of collisions can be ignored at the background level) and to a first-order perturbation equation 



dt 



dt 



+ 



a dx^ 



Hq 



djSf) 
dq 



~'~ dq 



a dx^ dt 



C[f]. 



(24) 



From Eqs. (fTTT) and (|23] ) we obtain that the temperature of the background photon temperature 



falls as 



T oc 1/a . 



From Eqs. (fT4]) and (1241) we obtain 



de q^de_ q^d^ ^-r\p,] 

dt a a dx'- dt 



(25) 



(26) 



the brightness equation. 



6 Thomson Scattering 

Photons scatter from charged particles. At the time of interest, these are electrons, protons, and helium 
nuclei. Since the scattering cross section is inversely proportional to the square of the mass of the charged 
particle, we need to consider just the electrons. In the non-relativistic (kinetic energies much below the 
electron mass) limit scattering of photons on electrons is called Thomson scattering. The differential 
cross section is 

do (Tt 3 / a „N 

Si = 4^4 + 

where 

aT = 0= 6.65 X lO^^V^ . (28) 

3 mj. 

The collision term Eq. (I26l ) is proportional to the electron density ng, where only free electrons 
count, not those already bound in atoms. If the electron fluid were in rest {Vf. = 0), the effect of scattering 
would be to isotropize the photon distribution, i.e., to damp all its higher moments. However, in the 
perturbed universe, there is a perturbation in the electron fluid velocity, equal to the baryon velocity 
perturbation, Ve = Vh, whose effect is to drag the photon velocity perturbation towards it. We skip the 
derivation of the collision term and just give the final form of the brightness equation: 



dQ q' de q' d^ 9* _ 

dt a 5x* a 5x* dt ""^"^^ 



eo-@{q) + q-Vb + lq'q^Ql' 



(29) 



The effect of the two first terms on the right-hand side (RHS) is to damp all multipoles of @{q), except 
00. The effect of the third term is to instead force = 30i towards Vh. The last term is due to the 
angular dependence of Eq. (l27l) . which has a quadrupole shape, and has the effect that the quadrupole of 
the photon distribution is not damped as fast as the other multipoles. 

The differential cross section actually depends also on photon polarization (Eq. (|27] ) is averaged 
over the two polarization directions). The quadrupolar angular dependence of this has the effect that Q2 
acts as a source of CMB polarization. In this lecture, however, we discuss just the generation of the CMB 
temperature anisotropy, not its polarization. 
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7 Line-of- Sight Integration 

The recombination of hydrogen had a dramatic effect on the brightness equation, since the density of 
free electrons rig dropped by many orders of magnitude. 

Before recombination, was large, forcing the term in the brackets in Eq. ( |29l ) to be very small. 
We can then make the tight-coupling approximation: 

e{q) = @Q + q-Vf, =^ v^ = 3ei= Vb] @^ = (30) 

After recombination. Tie was so small that most CMB photons have never scattered after recombi- 
nation. We can then make the collisionless approximation, and use the collisionless brightness equation 

de^^de_^q^d^_d^_^ ^^^^ 
dt a dx^ a dx'- dt 

The total derivative along a photon path can be written 



d _ d ^ d ^^^^ 
dt dt a dx^ 

and thus we get from (|3TI ) that along a photon path 

The quantity G + <I> is called the effective temperature perturbation, since it adds to the local temperature 
perturbation the effect of the gravitational red/blueshift from the local gravitational potential. 

Although in reality hydrogen recombination lasted tens of thousands of years, and an exact cal- 
culation has to follow this, we can get a good qualitative understanding of the CMB anisotropy by 
making the instantaneous decoupling approximation: we assume that recombination took place at t = 
(« 380 000 yr), and use the tight-coupling approximation for t < and the collisionless approximation 
for t > t*. 

We can then integrate Eq. (|33] ) along the photon path (line-of-sight integration), starting at time 
from the location where the photon last scattered, xis, to the present time to and the location Xohs where 
the photon is observed today: 

G>{to,Xob.,q) + ^{to,Xohs) = ie + '^){t^,xis,q) + \~dt ^ ~dt ) ^''^ ^^^^ 

= eo{t„xis)+^{t^,xis) + q-Vby + j (-^ + ^)^* (35) 

Apply now this result to a fixed observer, looking at all directions. The observed perturbation 
in the CMB temperature in direction h = {9,(j)) is given by Q{tQ,Xohs, —n), since the observer is 
looking against the photon direction. The term <I>(to, 2fobs)> which is just the gravitational potential of 
the observing site, does not depend on the direction looked at, and thus appears just an overall shift in 
the mean CMB temperature. This effect is smaller than the accuracy the mean CMB temperature has 
been measured with, and we ignore it. On the RHS we see four different contributions to the CMB 
temperature anisotropy: 

- The original temperature perturbation at the last scattering sphere, 0o(t*, x^s) = 

- The gravitational potential from which the CMB photons have to climb (or fall) from, <I>(t*, xj^), 
causing a gravitational red/blueshift of the radiation temperature 

- A Doppler effect q ■ Vhj = — n • Vby coming from the motion of the primordial baryon-photon fluid 
at the last scattering sphere 
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- An effect that comes from the time dependence of the metric perturbations along the photon path. 
If the gravitational potential $ does not depend on time, the redshift due to falling in it is canceled 
by the blueshift due to climbing from it. Thus gravitational potential along the photon path has a 
net effect only if it is time dependent, so that this cancellation is not exact. The same applies to the 
perturbation in the expansion rate. This effect is called the integrated Sachs-Wolfe (ISW) effect. 

The metric perturbations <I>, ^ are affected by all energy components, pi,, Pcdm> Pj^ Pv^ and there- 
fore we need the evolution equations for all of them, in addition to the Einstein equations from general 
relativity for the evolution of <I>, and ^. To obtain the quantities needed in Eq. (I35l) . we need to integrate 
these evolution equations starting from initial conditions specified at some time well before recombi- 
nation, when the universe was still radiation dominated, and all scales of interest were "outside the 
horizon", meaning that the Hubble distance was then smaller than these scales. 

These initial conditions for the perturbations are called primordial perturbations, and to obtain a 
theoretical prediction for the observed CMB anisotropy (as well as for the observed matter distribution 
today), we need a theory for the production of primordial perturbations. 



8 Primordial Perturbations 

The primordial perturbations were apparently produced by some random process. Therefore we only 
expect to predict their statistical properties. The current favorite scenario for their production is called 
inflation. Inflation refers to an accelerating expansion of the universe by a very large factor at very 
early times. Because the expansion is accelerating, the scale factor a grows much faster than the Hubble 
distance, causing perturbations to exit the horizon (their distance scale becomes larger than the Hubble 
length). During inflation microscopic scales were expanded to astronomical scales, and the primordial 
perturbations are produced from quantum fluctuations at these microscopic scales. 

There are many proposed theories where inflation can be realized, but in the simpler ones there is 
only one dynamically important independent quantity at that time, a scalar field if, called the inflaton. 
During inflation, the homogeneous background value ip{t) "rolls" slowly towards the minimum of the 
inflaton potential V{ip). All particles in the later universe are produced after inflation from the energy that 
was stored in the inflaton field during inflation, in a process called reheating. All perturbations arise from 
the inflaton perturbations 6ip. Since perturbations in all quantities originate from a single perturbation 
quantity, they are related to each other in a simple manner, i.e., the resulting primordial perturbations are 
adiabatic. 

This means, e.g., that the perturbations in the number densities of all particle species are the same 

J rii rij 4/9^ 

(The photon number density is related to temperature by a and the photon energy density by 
oc T^.) For baryons and CDM, pi = mini, so that 

5 Pi 6ni 3. 

Oi = — = = j0.y = dm- (37) 

Pi rii 

After inflation, all cosmological distances are much larger than the Hubble distance, and therefore 
perturbations at these scales do not have any dynamical evolution. These "superhorizon" perturbations 
are naturally described in terms of the associated spacetime curvature perturbation. We can define a "co- 
moving curvature perturbation" TZ{t, x) (related to ^> and ^) that stays constant in time at superhorizon 
scales (for adiabatic perturbations). 

The quantum fluctuations during inflation are a random process, and therefore we can not pre- 
dict individual values of Tl{x) from an inflation theory, but we can predict expectation values of the 
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magnitudes of the perturbations at different distance scales. These are given by the power spectrum 



rn{k) = X (|7^g|2) (38) 

where (•) denotes expectation value, and IZj: is the Fourier amplitude of TZ corresponding to wave vector 
k. Here V is the reference volume used to Fourier expand TZ{x), and its choice does not affect the result. 
The inflation prediction for V-jiik) is 

where Mpi is the Planck mass, is the value of the inflaton field when scale k exited the horizon 
{k = aH) during inflation, and 

e = «1. (40) 

The primordial perturbations produced in inflation are close to scale invariant, 

V-jiik) ~ const. , (41) 
since during inflation 99 and H change slowly, while the scale factor a grows rapidly. More accurately 

Vn{k) ~ Alk"""-^ where - 1 = -6e + 2r? (42) 

where is the spectral index of the perturbations and 

V" 

V = M].i—, |r/|<l. (43) 

For historical reasons, there is this —1 in the common definition of the spectral index n^, but the relevant 
quantity is 

dlnVn 

"'^ dink ' ^ ' 

which gives the scale dependence of the primordial perturbations. The slow-roll parameters e and ■>] 
depend on the inflation model, but they are always small in successful models. 

After inflation, as the universe gets older, the Hubble distance H^^ grows, faster than the scale 
factor a, and encompasses larger scales. At the photon decoupling time t*, the Hubble distance was 
^ 200 Mpc, corresponding to about 1° on the CMB sky. Thus at angles S> 1° we see superhorizon 
perturbations that are still in their primordial state, i.e., they have not evolved since they were produced. 

For these large scales, it is easy to derive an approximate prediction for the CMB anisotropy using 
Eqs. (l35l) and (l39l ). At the universe was already matter dominated by a factor of a few over radiation. 
Thus we approximate the total density perturbation 6 = 6p/p hy the matter density perturbation, 

6^6m = lS^ = 390 (45) 

It turns out that $ and ^ are constant in time in a matter-dominated universe. Thus the ISW effect gets 
a contribution only from early times after photon decoupling, when the universe is not yet completely 
matter dominated, and at late times when the universe became dominated by dark energy. Therefore the 
last term of Eq. (l35l) is subdominant and we ignore it in our approximation. For adiabatic perturbations, 
velocity perturbations are negligible at superhorizon scales, and thus we approximate Eq. (I35l) by 

^-^{e,(t>)^lS + ^ (46) 
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where the RHS refers to conditions on the last scattering sphere. 

Now we still need to relate 6, <I>, and TZ. In Newtonian gravity 



(47) 



Now our gravitational potential is due to density perturbations, so we replace p hy 6p = p6. Since we 
are using comoving coordinates, we replace by (l/a^)V^. Using Eq. dD we would then get 



2 f k 

3 \aH 



(48) 



in Fourier space. The correct result derived from general relativity is 

2' 



2 + 



k 

3 Va^ 



k ' 



(49) 



showing that we get the Newtonian result for subhorizon scales {k <^ aH). However, we are now 
discussing superhorizon scales, so we instead have the approximate result 



-2$ 



(50) 



From general relativity, $ and TZ are related by 

$ = - 

in a matter-dominated universe. 

Thus the CMB anisotropy at large scales is 

5T 



5'^ 



(51) 



+ $ = 



5 • 



(52) 



Before converting this result and Eq. ( [39l ) into a prediction of the statistical properties of CMB anisotropy, 
we need to discuss how the latter are described. 



9 CMB Angular Power Spectrum 

The observed CMB temperature variations form a function on a sphere (the celestial sphere). In general 
this refers just to the unit sphere of observation directions h, but in the approximation (l46l ) this corre- 
sponds to the last scattering sphere. A standard way to analyze functions on a sphere is the expansion in 
terms of spherical harmonics 

— iOA) = Y,aimYi„,{0,(l)), £ = 0,1,..., m = -i,-£ + l,...,£, (53) 

where the harmonic coefficients aim are obtained by 

aen. = I dQYl^ie, (/>)^(^, </>) . (54) 

This is analogous to the Fourier expansion of functions of three-dimensional space. The different multi- 
pole numbers i correspond to different angular scales, with a rough correspondence 



11 



The different m for a given i correspond to different patterns or orientations with the same angular scale. 

The aim depend linearly (through the linear physics of first order perturbation theory) on primor- 
dial perturbations. Since different Fourier modes evolve independently, their amplitudes are uncorrelated, 
and this lack of correlation is inherited by the multipole coefficients 

(a^mO^'m') = for £ ^ £' or m ^ m' , (56) 

The evolution of the perturbations is different for different distance scales, but the physics does not 
differentiate between directions, and therefore Eq. (1381 ) depends just on the magnitude k of k. The 
analogous property of the multipole coefficients is that the expectation values of their amplitudes depend 
on i only, not on m. This dependence is called the angular power spectrum: 

|2\ 



Ci = {\ae,n 



(57) 



From this one obtains that the expectation value of the square of the temperature perturbation is 
given by a sum over the angular power spectrum, 

2 



5T_ 
~T 



2 i + l 
47r 



(58) 



Fig. |4] shows the observed CMB angular power spectrum based on 7 years of measurements by 
the WMAP satellite |j9l. The strong peak near i ~ 200 corresponds to structure at 1° scale, which is 
prominent in Fig. ??. 




100 500 

Multipole moment / 

Fig. 4: The CMB angular power spectrum from 7 years of measurements by the WMAP satellite. Black dots with 
error bars represent the measurements and the red curve is a theoretical prediction from a best-fit cosmological 
model. The blue band represents the expected scatter (cosmic variance) due to the randomness of the generation 
of primordial perturbations. FromRefs. [6 ,7 |. Credit: WMAP Science Team. 



10 Large Scales 

We have a prediction for the three-dimensional power spectrum of primordial curvature perturbations, 
Eq. (l39l) . from inflation. On the other hand, Eq. (l52l) relates CMB anisotropy to the values of 

7^(fz,) = ^7^,v^^•^^'« (59) 

k 

on the last scattering sphere. To get from an expansion in terms of plane waves to a spherical harmonic 
expansion, we need the relation between 3-dimensional plane waves and spherical harmonics 

gifc.x,. = 47r ^ i'Mkxis)Ye,^> {x)Y^1^, {k) ■ (60) 

I'm' 
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Here xis is the coordinate distance to the last scattering sphere, and ji are spherical Bessel functions. 




Fig. |5] illustrates how a Fourier mode (a plane wave) contributes to different angular scales at 
different parts of the last scattering sphere. Thus a given wave number k contributes to many multipoles 
£, as given by je{kxig), but the maximum contribution is around 

i ~ kxis . (61) 

Now it is straightforward to calculate 

= ^ = ■■■ = % Y.^\n^\")^,{kxf = I / ^VnikMkx)' . (62) 

^ k 

For the case of a scale-invariant power spectrum, V-jz = const. , we can use the property 

dk / , , o 1 _ 



of spherical Bessel functions to arrive at 



^' = l5-WTT)^ ^^^^ 



or, using Eq. (I39l ). 



27r ' 25 6007r2M|; e ' ^ ' 

This result is the reason why the CMB angular power spectrum is customarily plotted as {£{£+1) / (27r))C£. 

Eq. (|65] ) should apply for large scales, ^ 1°, or £ <C 180. From Fig.|4]we read that the observed 
value is about lOOO/iK^/TQ « 1.3 x 10^^'^. This gives a constraint for inflation models 

pa 8 X 10"^Mi ^ (0.03Mpz)^ (66) 

where ip^ refers to the value of the inflaton field when the observed cosmological scales exited the 
horizon. Since e ^ 1, this gives an upper limit to the inflation energy scale 

V{^x)^^^ < OmUpi = 7 X 10^^ GeV . (67) 
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11 Smaller Scales 

For smaller scales, 6* < 1°, or ^ > 180, the perturbations enter the horizon before t*, and therefore they 
have had time to evolve from their primordial state before we observe them. We give only a qualitative 
discussion of the main features of the relevant physics. The gravitational attraction of the overdense 
regions causes the tightly coupled photon-baryon fluid to fall into their gravitational wells. However, this 
increases the radiation pressure, which eventually pushes it out. Thus the fluid begins to oscillate, moving 
in and out of the gravitational wells. The gravitational potential is dominated by cold dark matter, which 
does not feel the radiation pressure, and therefore does not participate in these acoustic oscillations. 
Different Fourier modes of the perturbations oscillate with different frequencies, the relation between 
the wave number k and (angular) frequency uj = Cgk given by the sound speed Cg in the baryon-photon 
fluid, 

cl = \-^ where = (68) 

The oscillation in the density of the baryon-photon fluid (or the temperature perturbation 9o pro- 
portional to it), is not symmetric, since it is biased by the gravitational potential of the CDM: the maxi- 
mum over-densities in the potential wells are always larger than the under-densities at the opposite phase 
of the oscillation. A proper calculation gives that the quantity that oscillates around zero is proportional 
to 9o + (1 + ii)<I>, so that we have 

Q^^ + {1 + R)<l>^(x cos cskt (69) 

if we ignore the expansion of the universe. The expansion causes R and Cs to change with time, so that 
we have to replace Cst in Eq. (l69l ) by the sound horizon 

rsit)^ t^^dt', (70) 



JO a{t') 

the coordinate distance travelled at sound speed by time t. 
The first two terms in Eq. (|35] ) are thus given by 

{Go + ^)f:{U) ^ -R^^iQ + A^coskrsiU) . (71) 

The amplitude Aj: is complicated to derive from TZj:, since it is affected by physics near the time of 
horizon entry, where neither superhorizon, nor subhorizon approximations apply; but for smaller scales 
it gets a notable boost by a gravitational driving effect due to the evolution of the gravitational potential 
<I> when the universe was not yet matter dominated: The baryon-photon fluid falls into an evolving 
gravitational well that is becoming weaker. Thus when it is kicked out by the radiation over-pressure, it 
flies out further than from where it came, boosting the oscillation amplitude. 

The most conspicuous feature of Eq. dTTI ) is its oscillatory dependence on scale k: It is maximal 
at those scales k, given by 

krg = rrnr , m = 1, 2, . . . (72) 

where the oscillations were just at their extrema when photons decoupled. Thus we see a strong CMB 
anisotropy at the corresponding angular scales 

" or £ ~ kxig = mvr — = mi a (73) 



kxig i^s{t*) 

where 

^ . = -77, 

is called the acoustic scale in multipole space and 



Ia = vr4^f^ EE ^ (74) 



(75) 
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is the sound horizon angle. 

This phenomenon is responsible for the oscillatory behavior in the Ci seen in Fig.|4] We get from 
the separation of the peaks the observed value for £a ~ 300, which is a tight constraint on cosmological 
parameters which (t* ) and xis depend on. 

12 Diffusion Damping 
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Fig. 6: The effective temperature, Go + contribution to the angular power spectrum C^, calculated both with 
and without the effect of diffusion damping. The spectrum is given for four different values of total matter density 
w„i, with baryonic matter density uj\, = 0.01. Figure and calculation by R. Keskitalo. From Ref. |i8J. 



The most important effect that we neglected in making the instantaneous decoupling approxima- 
tion is photon diffusion. During recombination the photon mean free path grows rapidly. While the 
photons are still scattering, the photons carry energy from one part of the fluid to another, and this effect 
acts towards homogenizing the fluid temperature over a distance scale related to the photon mean free 
path. This damps the temperature perturbations at the smaller scales. The effect on Ci is quite dramatic 
as can be seen in Fig. (6] 

13 Putting It All Together 

The effective temperature perturbation Bq + ^* is the dominant contribution to Q, but the other two 
terms in Eq. (|35] ) are important also. In acoustic oscillation the fluid velocity v^^^ oscillates too, but in a 
different phase: for a given Fourier mode, when the density is at the extremum, the fluid is momentarily 
at rest. Thus the contribution from the Doppler effect —n ■ v^y to Ce is also peaked, but at different i, 
and therefore it acts to partially fill the troughs between the peaks coming from ©o + ^. 

The is quadratic in 5T/T and therefore it also picks contributions from products of the three 
separate contributions ©o + ^, —n ■ Vh^, and / {d^ /dt + / dt). Since the first two oscillate in different 
phases, their cross-contribution is small, but the ISW effect, which alone is the smallest of the three, is 
correlated with the effective temperature perturbation, and therefore actually contributes more through 
its cross-term with © + see Fig. |7] 
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Fig. 7: The full Ci spectrum calculated for the cosmological model = 0.2, ojf, = 0.03, no dark energy, 
and Us — 1, and the different contributions to it. (The calculation involves some approximations which allow 
the description of Ce as just a sum of these contributions and is not accurate enough for estimating cosmological 
parameters from data.) Here Oi denotes the Doppler effect. Figure and calculation by R. Keskitalo. From Ref. ID. 



14 Cosmological Parameters 

We can compare the observed Ci to those predicted by different cosmological models to determine the 
values of the free parameters in those models. Other cosmological data can usefully supplement CMB 
observations, but accurate CMB observations are vital for this task. 

The simplest cosmological model that fits current data is the ACDM model, which has: 

- a flat background universe 

- primordial density perturbations with a constant spectral index, 

- which are adiabatic 

- no primordial gravitational waves 

- dark energy in the form of a cosmological constant A (vacuum energy with a constant density) 

- negligible (^ 1 eV) neutrino masses 

This model has 6 parameters: the amplitude Ag and spectral index Ug of primordial perturbations 
(see Eq. l42l) : the background densities of baryonic matter ujh, cold dark matter ujcdm, and the density 
of vacuum energy Q^; and the optical depth of the universe after recombination r. (In Figs. [6]and|7J 

= l^cdm + Wfe.) 

Here the densities are given in terms of density parameters 




(76) 



Pc 



where 



SttG 



1.88 X lO^^^/i^kg/m^ 



Pc = 




(77) 



16 



is the critical density for the universe required to make the background universe flat, and h = Hq/ (lOOkm/s/Mpc). 
For baryonic and cold dark matter we have further defined coi = Q.ih?, which give the average density of 
these components in the universe today as 

PiQ = u;iX 1.88 X lO^^^kg/m^ . (78) 

The optical depth gives the expectation number of scatterings per CMB photon after recombination. 
These are mainly due to free electrons that were liberated when the radiation from the first stars reionized 
the interstellar gas. 

15 WMAP Results 

The WMAP satelUte Q began observing in August 2001, and ceased operations nine years later. WMAP 
measured the microwave sky at 5 frequency bands centered from 23 to 94 GHz ITOl . To improve their 
sensitivity, the WMAP instruments were passively cooled with a solar shade and radiators to an operating 
temperature of 90 K. 

The WMAP team has so far published results based on the first 7 years of data. According to 
them 1 1 1], the values of the ACDM model parameters are 

As = 4.94 ± 0.05 X 10"^ 

Hs-l = -0.037 ±0.012 

UJb = 0.0226 ± 0.0005 

Wcdm = 0.112 ±0.004 

= 0.728 ±0.015 

T = 0.087 ±0.014 (79) 

Except for the spectral index — 1 and optical depth r, these are already quite accurate. However, we 
get this accuracy only by assuming the ACDM model. If we relax these assumptions, we can use the 
data to derive upper limits on deviations from the ACDM model, but at the same time the error bars for 
the ACDM parameters become larger, in some cases a lot. 

There are also a number of anomalies in the WMAP data that cast some doubt on whether our 
current understanding of the universe is correct lfT2l . These are features that appear statistically unlikely 
in the scenario for the generation of primordial perturbations we have described. These anomalies include 

- a low quadrupole: the observed C2 is well below the expectation value of from the best-fit model 
(see Fig. 

- "axis of evil": the pattern of the quadrupole and octupole parts of the CMB anisotropy are curiously 
aligned, having a common preferred direction [13] 

- cold spot: the cold region near the lower right edge in Fig. |2]is unusually deep for such a small 
feature El 

- north-south asymmetry: if one divides the celestial sphere along the ecliptic, the northern hemi- 
sphere (upper left in Fig. ^ has much less large-scale anisotropy than the southern hemisphere 
(lower right) 1 15 1 

16 Planck Mission 

The next major improvement in our knowledge of the CMB anisotropy and polarization will come from 
the Planck satellite LLiMJ, see Fig.[8l which is a mission of the European Space Agency (ESA). 
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Fig. 8: An artist's view of Planck. Between the primary and secondary mirrors is the focal plane unit (see Fig.fTOll. 
to where the mirrors direct the microwave radiation. These are protected from straylight by a baffle. Below and to 
the left of the baffle we see three layers of thermal shields that provide passive cooling. On the other side of the 
shields there is the warm service module. The solar panels and the antenna for communication with ESA's ground 
stations are on its other side and are not visible in this view. From Ref. Ojj. Credit: ESA/AOES Medialab. 



Planck was launched into space, together with the Herschel infrared space telescope, on an Ariane 
5 rocket from Europe's Spaceport in Kourou, French Guiana, on May 14th, 2009. Planck reached its 
orbit, around the 2nd Lagrange point (L2) of the Earth-Sun system, on July 2nd. See Fig.|9l 

Lagrange points are special points in celestial mechanics, where objects in free fall stay at the 
same relative position with respect to Earth and Sun. L2 is 1.5 million km from Earth in the anti-Sun 
direction. L2 itself is in Earth's shadow, and since Planck draws its electric power from solar panels, it 
must not go too near L2. Planck's orbit around L2 has a radius of 400 000 km. (WMAP is also in orbit 
around L2). 

Compared to WMAP, Planck is an improvement in three respects: 

- Planck has better angular resolution, 5 arcmin compared to WMAP's 14 arcmin 

- it has a wider frequency coverage, which is important for separating the CMB from foreground 
radiation 

- it has a higher sensitivity thanks to a much lower operating temperature 

The higher angular resolution makes it possible to measure to much higher multipoles. From 
Eq. (|55] ). while WMAP reached to £ ~ 800, Planck reaches to ^ ~ 2200 (actually one gets informa- 
tion also from somewhat higher multipoles, but with less accuracy). This is especially important for 
measuring the spectral index to higher accuracy. 

To cover the wider range of frequencies, Planck carries to separate instruments with two different 
detection technologies: The Low-Frequency Instrument (LFI) has radiometers for measurements at three 
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Fig. 9: Planck's cruise to and orbit around L2. From Ref. [1]. Credit: ESA - C. Carreau. 



frequency bands, centered at 30, 44, and 70 GHz. The High-Frequency Instrument (HFI) has bolometers 
for measurements at six frequency bands, centered at 100, 143, 217, 353, 545, and 857 GHz. See Fig.fTOl 

The higher frequencies are important for separating out the microwave radiation from interstellar 
dust, whose radiation increases with frequency. They are also important for measuring the Sunyaev- 
Zeldovich (SZ) effect, which is the upscattering of CMB photons to higher frequencies by the energetic 
electrons in the intergalactic space of galaxy clusters. In the direction of galaxy clusters, the SZ effect 
lowers the CMB intensity at frequencies below 217 GHz and raises it for frequencies above 217 GHz. 
See Fig. [TT] Although Planck is primarily a CMB mission, the full-sky observations at these many 
frequencies are also important for many astrophysical studies. 

To achieve the high sensitivity of its instruments, Planck carries a three-stage active cooling sys- 
tem. The first stage is a hydrogen sorption cooler to cool the LFI to 20 K. The second stage, the 4K 
cooler, is based on Joule-Thomson expansion of helium. In addition to being one stage of cooling the 
HFI, the 4K cooler also provides a reference load at 4 K temperature for the LFI radiometers. The fi- 
nal stage is a dilution cooler, which operates with helium-4 and helium-3, and is used to cool the HFI 
bolometers to 0. 1 K. Planck carries large tanks of helium-3 and helium-4, since after dilution these gases 
are vented to space. This limits the operating time of HFI, since this gas supply will eventually run out. 

The high sensitivity is especially important for polarization measurements, since the CMB polar- 
ization is at least an order of magnitude weaker than the temperature anisotropy. Except for the two 
highest frequencies, all Planck channels measure the polarization also. 

The noise of the LFI radiometers is dominated by low noise frequencies. To remove this low- 
frequency noise the radiometers continuously observe the 4K reference load together with the sky. The 
signal from these two sources are switched at 8192 Hz frequency between the diodes of the radiometers, 
and are afterwards separated again into a sky signal and a reference signal. In this way both signal 
streams have come through the same electronics, and contain the same low-frequency noise. Taking the 
difference between the two signals removes most of the it, leaving the sky signal and residual noise that 
is almost white (uncorrelated). 

Planck rotates at 1 rpm. The instruments are pointed at a direction about 85° away from the spin 
axis, scanning almost great circles on the sky. The spin axis is repointed by 2 arcmin about once per hour 
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Fig. 10: The feedhorns of the Planck receivers at the focal plane of the Planck telescope. The smaller horns at the 
center belong to the HFI. They are surrounded by the six 70 GHz, three 44 GHz, and two 30 GHz LFI feedhorns. 
FromRef. O]. Credit: ESA/AOES Medialab. 




70 6Hz 100 6Hi »aCHi 217 GHi 35SGHi B4SGHi 



Fig. 11: Abell 2319, a nearby cluster of galaxies, seen by seven different Planck channels. For frequencies below 
217 GHz, the cluster appears as a cold spot, for frequencies above 217 GHz as a hot spot. From Ref. [17] . Credit: 
ESA/LFI & HFI Consortia. 



to keep it pointed close to the anti-Sun direction. To observe also the regions near the ecliptic poles the 
repointing scheme actually follows a cycloid around the anti-Sun direction so that the spin axis always 
points 7.5° away from it. This way the whole sky is covered twice in one year. Repeated measurements 
of the same sky points at different times is used at the map-making stage 1.18.1 of the data analysis to 
remove residual correlated noise. 

Planck reached the 0.1 K temperature about 50 days after launch, making the HFI bolometers 
the coldest known place in space! After that followed a period of tuning of the instruments and perfor- 
mance verification. On August 12th, science observations began with a two-week "first-light survey" 
(see Fig. [Jl. Without break the observations continued into the originally planned 15-month "nominal 
survey". The whole sky was observed by June 2010. See Fig. [121 This nominal Planck mission ended 
on November 26th, 2010. ESA has, however, extended Planck operations, both with LFI and HFI, by 12 
months, until near the end of 2011. LFI does not need the dilution cooler, and in November 2010 ESA 
extended the Planck observation program yet by another year, for using LFI only. 
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The Planck one-year all-sky suruey ^GSa [ciESA,Hn=nnFi consortia, jui^joio 



Fig. 12: A full-sky image of the microwave sky based on the first year of Planck observations. The colors in the 
image are related to the dominating microwave frequencies at the different parts of the sky. This image is made 
combining different Planck frequencies so that the foreground-dominated ones have been given the most weight, 
and therefore the CMB, shown in red and orange is only visible far from the galactic plane. In other parts of the 
sky the image mainly shows emission from the gas and dust in our own galaxy. From Ref. G]. Credit: ESA/ LFI 
& HFI Consortia. 



ESA has granted the Planck Collaboration a two-year proprietary period for data analysis and 
deriving science results, after which the data will be made public. Thus the data from the nominal 15- 
month mission will be released near end of 2012, and the main cosmology results are expected to be 
published then also. Some early results on the astrophysics of foregrounds will already be published in 
January 2011, including the Early Release Compact Source Catalog. 

17 Cosmology from Planck 

Planck will extract almost all the available information from CMB temperature anisotropy, limited by 
cosmic variance and the ability to separate foregrounds from the CMB. Planck will also provide an ac- 
curate spectrum of CMB polarization for the first time, although Planck is not optimized for polarization 
measurements, leaving room for a future CMB space mission focusing on polarization. 

The higher sensitivity and resolution of Planck will lead to a major improvement in the accuracy 
of the determination of cosmological parameters. The main cosmological interest is in the nature of 
primordial perturbations. This is the key to the mechanism for the origin of structure in the universe, 
which is responsible for the existence of galaxies, stars, and planets. This mechanism is related to very 
high-energy physics whose study is beyond the reach of Earth-based accelerators. 

A better determination of the spectral index Ug of primordial perturbations is already important for 
selecting among candidate inflation theories. But even more important is the ability of Planck to probe 
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Fig. 13: An illustration of the expected improvement in cosmological parameter determination for a 7-parameter 
model: ACDM supplemented with an additional parameter rimn = dug/d In k, allowing the spectral index to vary 
with distance scale. Blue contours show forecasts for WMAP after 4 years of observation and red contours show 
results for Planck after 1 year of observations. From Ref. IIT9l . 



additional cosmological parameters beyond the simple ACDM model (see Fig. [T3] ). Many inflation 
models predict also the production of gravitational waves, also called tensor perturbations, in addition 
to the density (scalar) perturbations. The upper limit from WMAP on the ratio of tensor perturbation 
amplitude to scalar perturbation amplitude is r < 0.24 (95% CL) [11], assuming that this is the only 
extension to the ACDM model. The effect of tensor perturbations on CMB temperature anisotropy is 
somewhat similar to some other cosmological parameters, but it causes an unambiguous signal in CMB 
polarization. 

The polarization field on the celestial sphere can be divided into an E-mode (curl-free part) and 
a B-mode (source-free). To first order in perturbation theory, scalar perturbations produce only E-mode 
polarization. Therefore a detection of B-mode polarization at relatively large scales £ < 150 is clear 
evidence of primordial tensor perturbations. At smaller scales, second-order effects convert a part of 
E-mode polarization into B-mode. The WMAP upper limit to tensor perturbations comes from CMB 
temperature anisotropy; the B-mode is beyond the reach of WMAP polarization sensitivity. Planck is 
sensitive enough to detect B-mode polarization coming from r < 0.1 [19.| . 

More complicated inflation models have additional signatures that are not included in the ACDM 
model, and Planck will be looking for these. They include deviations from the adiabaticity of primordial 
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perturbations, and the deviation of their statistics from Gaussianity. 

Another important thing expected from Planck is a check on the large-scale anomalies in WMAP 
data (Sec. [BJl. 
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